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ABSTRACT 

Context. We report on a successful, simultaneous observation and modelling of the millimeter (mm) to near-infrared (NIR) flare 
emission of the Sgr A* counterpart associated with the supermassive (4x10^ M© ) black hole at the Galactic centre (GC). We present 
a mm/sub-mm light curve of Sgr A* with one of the highest quality continuous time coverages. 
Aims. We study and model the physical processes giving rise to the variable emission of Sgr A*. 

Methods. Our non-relativistic modelling is based on simultaneous observations carried out in May 2007 and 2008, using the NACO 
adaptive optics (AO) instrument at the ESO's VLT and the mm telescope arrays CARMA in California, ATCA in Australia, and the 
30 m IRAM telescope in Spain. We emphasize the importance of multi- wavelength simultaneous fitting as a tool for imposing adequate 
constraints on the flare modelling. We present a new method for obtaining concatenated light curves of the compact mm-source Sgr A* 
from single dish telescopes and interferometers in the presence of significant flux density contributions from an extended and only 
partially resolved source. 

Results. The observations detect flaring activity in both the mm domain and the NIR. Inspection and modelling of the light curves 
show that in the case of the flare event on 17 May 2007, the mm emission follows the NIR flare emission with a delay of 1.5+0.5 
hours. On 15 May 2007, the NIR flare emission is also followed by elevated mm-emission. We explain the flare emission delay by 
an adiabatic expansion of source components. For two other NIR flares, we can only provide an upper limit to any accompanying 
mm-emission of about 0.2 Jy. The derived physical quantities that describe the flare emission give a source component expansion 
speed of Vexp ~ 0.005c - 0.017c, source sizes of about one Schwarzschild radius, flux densities of a few Janskys, and spectral indices 
of Qf=0.6 to 1.3. These source components peak in the THz regime. 

Conclusions. These parameters suggest that either the adiabatically expanding source components have a bulk motion greater than 
Vexp or the expanding material contributes to a corona or disk, confined to the immediate surroundings of Sgr A*. Applying the flux 
density values or limits in the mm- and X-ray domain to the observed flare events constrains the turnover frequency of the synchrotron 
components that are on average not lower than about 1 THz, such that the optically thick peak flux densities at or below these turnover 
frequencies do not exceed, on average, about ~1 Jy. 

Key words, black hole physics, infrared: general, accretion, accretion disks, radio, Galaxy: centre, nucleus 



1. Introduction 

Sgr A*, the compact non-thermal radio and infrared source at 
the centre of the Milky Way galaxy (~8kpc away) is known 
to be associated with a supermassive black hole (SMBH) of 
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mass -4 x10^ Mq (lEcka rt & Genzel||l996[ [Oenzel et al]| 19971 
2000; G hez et al.||199 8, 2000, 2004b a, 2005; Eckart et al] 
2002 ; Schodel et al.||2002, ,2003, Eisenhauer et al.,20Q 3, 2005} 
Gille ssen et al.|2009]|Ghez et al.|2Q09| ). The close'proximity of 
Sgr A* makes it ideal for studying the evolution and physics of 
SMBHs located in the nuclei of galaxies. The SMBH radiates 
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far below its Eddington luminosity at all wavelengths, partly be- 
cause of its low observed accretion rate. Eor Sgr A*, we assume 
that Rs=2Rg=2GM/c^^9 yuas, Rg being one Schwarzschild radius 
and Rg the gravitational radius of the SMBH. 

Evidence of flaring activity occurring from a few hours to 
days has been found from vari ability studies ranging from the 
radio to sub-mm wavelengths (Bower et al. 2002 ^ 'Herrnstein 
let al |2004[ [Zhao e t al. 2003 , 2004 ; Mauerhan et al . 2005 ). There 
is also evidence that variations in radio/sub-nmi emission are 
linked to NIR/X-ray flares, with the radio/sub-nmi flares occur- 
ring after a delay of ~100 minutes after the NIR/X-ray flares 
dEckart et al.||2008bl[2006a|[2Q04t|Marrone et al.|[20Q81|Yusef-| 
IZadeh et al.|2 008|). 

These flares have been explained with a synchrotron self 
Compton (SSC) model that involves up-scattered sub-mm pho- 
tons from a compact (inferred from short flare timescales) source 
component (e.g., Eckart et al. 2004, Eckart et al. 2006a). The X- 
ray emission is caused by the inverse Compton scattering of the 
THz-peaked flare spectrum by relativistic electrons. Both syn- 
chrotron and SSC mechanisms contribute to NIR flux density. 
Adiabatic expansion of the source components then gives rise to 
flares in the mm/sub-mm regimes ([Eckart et aL||2006a[ [Yusef-[ 



[Zadeh et al.|2006a|b[ [20081 [Marrone et al.|2008D . 

Based on the assigned VLT time, we organized an extensive 
multi-frequency campaign in May 2007 and May 2008 which 
included millimeter to NIR observations at single telescopes and 
interferometers around the worl d. We present data from observa- 
tions of Sgr A* using CARMApTlBock et al.|2006 ), ATCA[5 and 
the MAMBO bolometer at the IRAlvj^SO m telescope in the mm 
regime and NIR data from the ESO VLT. We detected simulta- 
neous emission in the NIR and mm-regimes using the CARMA 
and ESO VLT telescopes. The main results are: 2 bright NIR 
flares (16mJy and lOmJy), a ~0.4Jymm flare, and a possible 
weaker third flare in the NIR and mm in the combined ATCA, 
CARMA, IRAM 30 m mm/sub-mm light curve from May 2007, 
and a bright NIR flare in May 2008 also covered by the CARMA 
3 mm observations. We present updated versions of light curves 
first presented in [Kunneriath et al^ ( [200 8), a more detailed de- 
scription of the methods used to obtain them, and physical mod- 
els to explain the flaring activity. The observations and data re- 
duction are described in Sect. 2, an outline of the flare modelling 
in Sect. 3 and a discussion and summary follow in Sects. 4 and 
5. 



surroundi ng the Galactic centr e at a distance of about 1.5- 4 pc 
(see, e.g.,'Guesten et al.[( pW7l ) or [Christopher et al.[ ( [20Q5] )). In 
May 2007 and 2008, global coordinated multiwavelength obser- 
vations were carried out in the NIR and mm regimes to study the 
variability of Sgr A*. 

2.1. The mm data 

We observed the GC at 100 and 86 GHz (3 and 3.5 mm wave- 
length) with the two mm-arrays CARMA and ATCA, respec- 
tively. In addition, we observed with the MAMBO 2 bolome- 
ter at the IRAM 30 m-telescope at a wavelength of 1.2 mm. 
CARMA is located in Cedar Flat, Eastern California, and con- 
sists of 15 antennas (6 x 10.4 m and 9 x 6.1m). The quasar 
source 3C273 was used for bandpass calibration while 1733- 
130 and Uranus were used for phase and amplitude calibration 
of Sgr A* data, respectively. The Australia Telescope Compact 
Array (ATCA), at the Paul Wild Observatory, is an array of six 
22-m telescopes located in Australia. Calibrator sources 1253- 
055, 1921-293, and Uranus were used for bandpass, phase, and 
amplitude calibration. The interferometer data were mapped us- 
ing the Miriad interferometric data reduction package. Details of 
the observation are given in Table [T] 

The Max-Planck Millimeter Bolometer (MAMBO 2) array 
is installed at the IRAM 30 m telescope on Pico Veleta, Spain. 
The 37 channel array of the precursor instrument MAMBO has 
been successfully used by many observers since the end of 1998. 
The bolometer data was reduced using the bolometer array data 
reduction, analysis, and handling software package, the BoA 
(Bolometer Data Analysis). 

To correct for extended flux contributions in the interferom- 
eter data, we extracted visibilities from two orthogonal pairs of 
the longest baselines for the two arrays, and subtracted for each 
baseline the median baseline and time dependent visibility trend 



Db{t) = //epoch(4(0) 



(1) 



from each visibility data set dh(t). Here the operator yUepoch repre- 
sents the median over all epochs. The time-dependent difl'eren- 
tial visibilities S (t) and their uncertainties SS (t) were then cal- 
culated via 



S(t)=jit(dt(t)-Dt(t)) 



(2) 



2. Observations and data reduction 

Interferometric observations in the mm/sub-mm wavelength do- 
main are especially well suited to difl'erentiating the flux den- 
sity contribution of SgrA* from the thermal emission of the 
circumnuclear disk (CND, a ring-like structure of gas and dust 

^ Support for CARMA construction was derived from the states 
of California, Illinois, and Maryland, the Gordon and Betty Moore 
Foundation, the Kenneth T. and Eileen L. Norris Foundation, the 
Associates of the California Institute of Technology, and the National 
Science Foundation. Ongoing CARMA development and operations 
are supported by the National Science Foundation under a cooperative 
agreement, and by the CARMA partner universities. 

^ ATCA is operated by the Australia Telescope National Facility, 
a division of CSIRO, which also includes the ATNF Headquarters 
at Marsfield in Sydney, the Parkes Observatory and the Mopra 
Observatory near Coonabarabran. 

^ The IRAM 30 m millimeter telescope is operated by the Institute 
for Radioastronomy at millimeter wavelengths - Granada, Spain, and 
Grenoble, France. 



and 

SS(t)=fit(S(t)-(dt(t)-Dtm . (3) 

Here the operator is the median over difl'erent baselines. 
Hence, SS(t) is the median of the deviation from the median 
flux S(t). The visibilities had been calibrated via intermittent 
flux reference observations. We attribute the residual flux den- 
sity dips/excesses to variations in the intrinsic flux density of 
Sgr A*. The combined light curve from all telescopes is shown 
in Fig. 1. Here we can combine the data from difl'erent frequen- 
cies based on the assumption that the spectral index of Sgr A* 
does not change significantly during the flux density variations 
between 86 and 250 GHz. In this case, the flux density variations 
are frequency independent. 

The influence of flux spectral index variations can be calcu- 
lated in the following way. If S i and S 2 are the flux densities at 
the frequencies vi and y2, then a change from spectral index a to 
a -\- Aa results in a new flux density value S [ and with 

log Si = log S 2-^ Q: (log Vi - log y2) (4) 
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lelescope 


Instrument/ Array 


A 


U 1 and ju 


U 1 and JU 


Observing ID 






Start Time 


Stop Time 


ATCA 


H214 


3.5 mm 


2007 15 May 07:36:05 
JD 2454235.81672 


15 May 22:57:02.5 
JD 2454236.45628 


IRAM 30m 


A T A A 11 A 

MAMBO bolometer 


1.2 mm 


2007 16 May 00:20:42 
JD 2454236.51437 


16 May 04:18:56 
JD 2454236.67981 


CARMA 


D array 


3.0 mm 


2007 16 May 07:43:31.3 
JD 2454236.82189 


1 A T 1 O T7 r\'~i o 

16 May 13:27:07.8 
JD 2454237.06051 


ATCA 


H214 


3.5 mm 


2007 16 May 09:31:57.5 
JD 2454236.89719 


16 May 22:21:22.5 
JD 2454237.43151 


IRAM 30m 


A T A A 11 A 

MAMBO bolometer 


1.2 mm 


2007 17 May 00:14:39 
JD 2454237.51017 


17 May 04:28:13 
JD 2454237.68626 


A T~> A T A 

CARMA 


D array 


3.0 mm 


2007 17 May 07:22:46 
JD 2454237.80748 


17 May 13:20:43 
JD 2454238.05605 


ATCA 


H214 


3.5 mm 


2007 17 May 09:47:17.5 
JD 2454237.90784 


17 May 18:22:32.5 
JD 2454238.26565 


IRAM 30 m 


A /r A A 11 A 

MAMBO bolometer 


1.2 mm 


2007 18 May 00:12:57 
JD 2454238.50899 


18 May 04:23:03 
JD 2454238.68267 


A T> A /r A 

CARMA 


D array 


3.0 mm 


^r\r\'~i 1 o A T r\'~i 0001 

2007 18 May 07:33:21.8 
JD 2454238.81484 


18 May 13:19:20.3 
JD 2454239.05510 


ATCA 


H214 


3.5 mm 


2007 18 May 10:11:27.5 
JD 2454238.92462 


18 May 22:23:47.5 
JD 2454236.67981 




D array 


3.0 mm 


9007 1Q A/Tq^/ 07-^0-^9 ^ 

JD 2454239.81311 


1 Q A/Tq^/ 1 1 ^-49 8 

i^* iviay ij.ij.^z.o 
JD 2454240.05258 


ATCA 


H214 


3.5 mm 


2007 19 May 09:30:57.5 
JD 2454239.89650 


19 May 22:13:07.5 
JD 2454240.42578 


CARMA 


C array 


3.0 mm 


2008 26 May 06:45:13 
JD 2454612.78140 


26 May 13:10:43 
JD 2454613.04911 



Table 1. Log of the mm- and sub-mm observations. 



Telescope 


Instrument/Array 




X 


UT and JD 




UT and JD 


Observing ID 








Start Time 




Stop Time 


VLT UT 4 


NACO 


2 


2 jim 


2007 15 May 05:29 
JD 2454235.72847 


00 


15 May 09:42:00 
JD 2454235.90417 


VLT UT 4 


NACO 


2 


2 jim 


2007 16 May 04:47 
JD 2454236.69956 


22 


16 May 07:54:41 
JD 2454236.82339 


VLT UT 4 


NACO 


2 


2 jim 


2007 17 May 04:42 
JD 2454235.69600 


14 


17 May 09:34:40 
JD 2454235.89907 


VLT UT 4 


NACO 


2 


2 fim 


2007 19 May 04:55 
JD 2454239.70486 


00 


19 May 09:28:22 
JD 2454239.89470 


VLT UT 4 


NACO 


3 


8 jim 


2007 15 May 10:05 
JD 2454235.92069 


48 


15 May 10:32:40 
JD 2454235.93935 


VLT UT 4 


NACO 


3 


8 jim 


2007 16 May 08:28 
JD 2454236.85317 


34 


16 May 10:44:04 
JD 2454236.94727 


VLT UT 4 


NACO 


3 


8 jim 


2007 17 May 10:16 
JD 2454237.92806 


24 


17 May 10:28:00 
JD 2454237.93611 


VLT UT 4 


NACO 


3 


8 yum 


2007 18 May 06:03 
JD 2454238.75238 


26 


18 May 10:29:05 
JD 2454238.93686 


VLT UT 4 


NACO 


3 


8 jim 


2007 22 May 04:57 
JD 2454242.70630 


04 


22 May 06:22:01 
JD 2454242.76529 


VLT UT 4 


NACO 


3 


8 //m 


2007 23 May 04:39 
JD 2454243.69384 


08 


23 May 10:36:24 
JD 2454243.94194 


VLT UT 4 


NACO 


3 


8 yL/m 


2008 26 May 05:42 
JD 2454612.73814 


55 


26 May 10:37:39 
JD 2454612.94281 



Table 2. Log of the near-infrared observations. 



log S[ = log ^2 + (o^ + Aa)(log vi - log 



we find 



V2 



or for 5; =/^i 



logf = Aalog— . 

V2 



(5) 



(6) 



(7) 



We estimate a spectral index variation of Aa = 0.07 between 
10^^ Hz and 10^^ Hz from the compilation of radio to sub-mm 
measurements by Mar koff^ et al.^ (2001) and Zhao et al.| p003| ). 
The maximum spectral range that we cover with CARMA, 
ATCA, and the IRAM 30 m is 86 GHz to 250 GHz. In total, 
this gives a maximum expected flux density variation of 8% 
in the light curve due to spectral index variations between the 
ATCA, CARMA, and the IRAM 30 m data. The bulk of the 
data in which we detected the main flux density excursions is 
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M,a^17 



May 19 



12 



12 



12 



12 



12 24 




VLBA 86 GHz 
VLBA 43 GHz 



NIR K & L band 



6 12 18 24 30 36 42 48 



54 60 66 72 78 84 90 96 102 108 114 120 126 132 138 144 
UT [hours] 



Fig. 1. Combined differential light curve of Sgr A* in the mm/sub-mm domain for the May 2007 observing run. The MAMBO 2 
bolometer at the IRAM 30 m- telescope was operated at a central wavelength of 1.2 mm (250 GHz). The CARMA data were centred 
on 100 GHz and the ATCA data on 86 GHz. We also show the daily 7 mm and 3 mm flux density averages from simultaneous VLBA 
observations ( |Lu et al.|2008[|20(j9] ). In the lower part of the plot, we show the combined NIR light curve consisting of K-band and 
L'-band light curves (see Appendix for individual light curves). The time axis is labelled with UT hours starting at 00 h on May 15. 



covered by the 100 GHz CARMA and the 86 GHz ATCA data. 
The CARMA and ATCA data follow the VLT data set, while the 
30 m data do not detect significant flux density excursions pre- 
ceding the VLT data. Here the maximum expected flux density 
variation in the light curve due to spectral index variations can 
only be 0.005 Jy, which is smaller than 6S(t). For the CARMA 
and IRAM 30 m data, the median flux density variation for each 
epoch (day) is SS (t) ~ 0. 1 Jy, and for the ATCA data we find that 
SS(t) ~ 0.2 Jy. These values represent the approximate range in 
which the datasets can be freely shifted in flux density. 

The light curve in Fig. 1 contains two peaks, on May 15 and 
17 (there is a weaker, third possible peak on May 19). In Fig. 
1, we also show the daily flux density averages of the 7 mm 



VLBA observations that were conducted in parallel ( |Lu et al. 
|2QQ8l [2009,. Kunneriath et al.,2 008). The VLBA data follow the 
overall trend of the combined CARMA/ATCA/30 m light-curve 
very well. The NIR coverage in the K and L'-band is given in 
the lower panel of Fig. [T] (see Appendix for a description of the 
scaling and the individual light curves for each day). 

To verify the mm/sub-mm flux density variations and in par- 
ticular the flare detected on 17 May by CARMA, we produced 
residual maps from the four individual tracks obtained with the 
CARMA array, since CARMA has the best combination of high 
resolution and signal-to-noise ratio of all the mm-telescopes we 
use. These maps shown in Fig. [2] were computed by subtract- 



ing the mean of all 4 maps from the maps of individual epochs. 
The rms noise in these difl'erential maps is of the order of 0.1 Jy 
per beam. This procedure clearly represents the trend shown in 
the combined light curve and detects the excess flux density of 
0.4 Jy detected on May 17 and a slightly positive flux density on 
May 19. On the 16th and 18th May 2007, a mixture of negative 
and positive flux densities in the difl'erence maps corresponds to 
the non detection of flaring activity in the difl'erential light curve 
shown in Fig.[T] 



2.2. The NIR data 

The NIR data were taken in the K^-band (2.18yum, FWHM 
0.35 fim) and in the L'-band (3.8 jam, FWHM 0.62 yum) at ESO's 
Very Large Telescope (VLT) with the NACO infrared camera 
and adaptive optics (AO) system ( Rousset et al.||2003{ |Lenzen 
|et al.|[2QQ3] ) in Chile. For the 3.8 //m observations, the inte- 
gration times were NDITxDIT=200x0.2=40 seconds (DIT is 
Detector Integration Time, and NDIT is the number of DITs). At 
2.2yum, the integration times were NDITxDIT=4x 10=40 sec- 
onds. Additional details of the observation are given in Table [2 
The diflraction limit in K-band is about 60 mas. The infrarec 
wavefront sensor of NAOS was used to lock the AO loop on 
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Fig. 2. Difference maps at 3 mm of a 4C diameter region (shown by the dashed white line) centred on Sgr A*, obtained from the 
difference between full synthesis maps of the individual days of CARMA observations and the full CARMA data set as described 
in the text. The maps have been plotted with the same color-coding table. The figure shows that the flux density variations that are 
evident from the differential light curves (see Fig.[T]) can also be seen in the maps constructed from the corresponding data. 



the NIR bright (K-band magnitude ~6.5) supergiant IRS 7, lo- 
cated about 5.6'' north of Sgr A*. Therefore the AO was able to 
provide a stable correction with a high Strehl ratio (of the order 
of 50%). After correcting the images for bad pixels, and sub- 
tracting the sky and flat fielding, the point spread function (PSF) 
was extracted from each individual image using StarFinder 
( Diolaiti et al.||2QQQ ) and a Lucy-Richardson (LR) deconvolu- 
tion was applied. A Gaussian beam with FWHM corresponding 
to the respective wavelength (-60 mas at 2.2yum and -104 mas 
at 3.8 yum) was used to restore the beam. To obtain a flux den- 
sity in Jansky, we perform aperture photometry of the sources 
with circular apertures of radius 52 mas, apply extinction cor- 
rections of A^=2.8 and A/^/ = 1.8, and calibrate the flux with flux 
densities of known sources, which resulted in K- and L'-band 
flux densities of S2 of 22+1 mJy and 9+1 mJy, respectively. The 
measurement uncertainties were obtained from the reference star 
S2. 



2.3. Flare events 

2.3.1 . The 15 May 2007 flare event 

The first NIR flare on May 15 detected during the multi wave- 
length campaign (see also Eckart et al. 2008a ) preceded the com- 
bined mm/sub-mm monitoring and the first maximum detected 
therein by nearly 3 hours. The May 15 2007 flare shows 2 po- 
larized bright sub-flares centred on about 07:29 UT and 07:51 
UT, respectively, after the start of observations (see Fig. [12] in 



Appendix). We define sub-flares to be the shorter flux excur- 
sions superimposed on the longer main underlying flare. The 
time difference between the two sub-flares is 22 minutes which 
is fully consistent with pre viously reported sub-flare separations. 
Zamaninasab et al. ( 2010| ) demonstrated that these highly polar- 
ized sub-flares are part of the flare structure that is significant 
compared to the randomly polarized red-noise. 

The second sub-flar e shows substructure that is interpreted 
by [Eckart et al. (2008a I in the framework of spot evolution. The 
flux density between the 2 sub-flares does not reach the emission 
level well before and after (i.e., <70 and >170 minutes into the 
observations) the flare, probably because of fore- or background 
stellar flux density contributions (star S 17 ). Additional details o f 
the May 15 2007 NIR flare were given by [Eckart et"aLl ( |2008a| ). 



For the flux density limit re ached in low lumin osity states of 
Sgr A*, the reader is referred to! Sabha et al.[ ( [201Q| ) and [Do et"aL 
(2009). The brief L'-band observations at 10 UT show a NIR 
sign al that is a factor of 2 times stronger than at other times (see 
Sect. [331] ). 



2.3.2. The 17 May 2007 flare event 



The NIR flare on 17 May had overlap with the CARMA obser- 
vations, which indicate that a 0.4 Jy 3 mm flare followed the NIR 
flare. 

Figure |3] gives the updated differential light curve obtained 
from the four longest baselines of CARMA and ATCA, along 
with the NIR light curve and a cross-correlation between the two 
showing a time lag of 1.5+0.5 hours, with a significant, broad 
and positive peak corresponding to the flares found at both wave- 
lengths. The presence of negative power in the cross-correlation 
is most likely due to a combination of two effects: (1) The signal- 
to-noise ratio of the differential ATCA data is lower than that 
of the CARMA data; and (2) the negative power may also be 
indicative of residual power on timescales of the length of the 
(sub-)n[mi-datasets. Some amount of negative power is expected 
since we are dealing with differential light curves that measure 
the variable flux with respect to a median flux level over the 
length of the corresponding dataset. A significant, broad, and 
positive cross-correlation peak centred on a time lag of -1.5 
hours is also present if the correlation is only performed with 
the K-band and CARMA data. Figure |4] shows that the calibra- 
tor and source fluxes are not correlated for the CARMA data. 
Here the cross-correlation was carried out only for the CARMA 
data since the ATCA array did not observe the same calibrator 
interleaved with the Sgr A* data. 

The NIR K-band data show variable emission, with four sub- 
flares starting from about 7:00 UT and lasting until the end of 
the observations at 9:12 UT. The mm-data started at about 8:00 
UT and had ended by 10:00 UT. In the mm domain, the variable 
emission is dominated by a single flare that is slightly wider than 
the individual NIR sub-flares. 
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Fig. 3. The 3 mm CARMA and ATCA (top red and blue, re- ^^^ .^^ CARMA SgrA*(top) and calibrator (1733-130) 

spectively) and K- and L'-band NIR (middle red and blue, re- (middle) light curves from 17 May 2007, with the cross- 

spectively) hght curves of the 17 May 2007 flare, with the cross- correlation (bottom) between the two showing zero correlation 

correlation (bottom) between the two showing that the mm flare between the calibrator and the source fluxes, 
follows the NIR flare with a time lag of 1.5+0.5 hours. ZDCF 
stands for Z-transformed discrete correlation function. 



2.3.3. The 19 May 2007 flare event 

The flux density excursions on 19 May are only partially con- 
strained by the CARMA millimeter measurements. Two flare 
events of about 3.5 mJy and 1.1 mJy peak flux density are sep- 
arated by about 3 hours. The gap in the observations at a time 
during which the first event reached its peak flux density (see 
Appendix for individual light curve) implies that the true peak 
was missed and that the flare was probably brighter than 3.5 mJy. 
There is an overlap between the NIR and the mm-measurements 
of about 1.5 hours with no detection of a flare event in the mm- 
wavelength domain brighter than 0.1-0.2 Jy. 



2.3.4. The 26 May 2008 flare event 

The 26 May 2008 NIR flare in the L'-band was covered by the 
CARMA observations, as shown in Fig. [5] The L'-band flare 
shows a flux density increase of about 70mJy, followed by a 
plateau with almost constant flux density at a 40 mJy level. After 
a total flare duration of 200 minutes, the flux density reaches 
the original level of 20 mJy. The CARMA data were taken un- 
der moderate to poor weather conditions. From 07:10 UT to 
11:10 UT, the source visibility was comparable to the levels 
reached in the 2007 observing session. No flare was detected 
with a flux density limit of about 0.2 Jy. Starting at 11:10 UT, 
the weather conditions deteriorated resulting in strong coherence 
losses. 
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Fig. 5. 3.8 yum NIR flare (top) covered by 3 mm CARMA data 
(bottom) on the 26 May, 2008. The time coverage of the L'-band 
flare is from 05:42:55 to 10:37:39 UT. 



3. Flare analysis 

The seven coordinated Sgr A* measurem ents reported so fa r that 
include sub-mm/mm data ([Eckart et~al . 2006a, 2008b ; Yusef- 
Zadeh et al.|[20Q6bt [Marrone et ari|2008, .Yusef-Zadeh et al. 
2009 including this work) have shown that the observed sub- 
millimeter/mm flares follow the largest event observed at shorter 
wavelengths (NIR/X-ray; see detailed discussion in Marrone 
et al.[ f2008)) with a delay of -1.5+0.5 hours. We therefore as- 
sume that the millimeter flares presented here are related to the 
observed IR flare events. This flaring activity can be explained 
using a multi-component relativistic spot/disk model and an adi- 
abatic expansion model in combination with a synchrotron self 
Compton (SSC) formalism. While this model has been presented 
before (see references above), it is essential to test new datasets 
against it to improve our understanding of the variable emission 
from Sgr A*. 



3.1. Relativistic spot/disk modelling of NIR flares 

Multiple hot spots revolving around a black hole in an accre- 
tion disk in AGNs or galactic black holes can give rise to light 
curves in the X-ray and NIR regimes whose power spectral den- 
sity (PSD) is described by a broken power law with a slope sim- 
ilar to that of red noise processes ( Armitage & Reyn olds 2003 ; 
Pechacek et al.|2008t|Do et al.|200'8, ). Simulations show that po- 
larized light curves exhibit behaviou r associated with tensing o f 
hot spots orbiting around a SMBH ( Zamaninasab et alT]|2010| ). 
Previously observed NIR light curves of Sgr A* have been suc- 
cessfully modelled by a multi-component hot spot model involv- 
ing source components orbiting around the SMBH in a tempo- 



This relativistic model consists of source components that 
are a mix of synchrotron and SSC components with an optically 
thin spectral index a and relativistic electrons with y^-lO^, using 
the formahsm des cribed by[Gould[(p ^979) and Mar scher[ ( |1983[ ). 
The KY-code by Dovci ak et al.[ (>2()()4) produces light amplifica- 
tion curves for each individual component orbiting the SMBH 
taking into account relativistic eff'ects, and by combining these 
curves with the SSC model, we can estima te the SSC X-ray and 
NIR flux densities and the magnetic field ( [Eckart e t al.||2008a|). 
Figure [6] presents simulations of the accretion disk with multi- 
ple spots revolving around the SMB H, and Fig.[7]gives the light 
curves produced by these simulations. Zamaninasab et al.[(2010|) 
describe the modelling and simulations shown in Figs. [6 and|7j 
in greater detail. 

These light curves may be the result of a single hot spot that 
dominates the disk for several orbital periods. If the spot were to 
sink (as may occur) towards the centre while orbiting, this would 
cause a quasi-periodic modulation of the light curve for the time 
the emission of the spot dominates the red-noise emission of the 
remaining disk. This modulation may also be produced by sev- 
eral diff'erent hot spots that have individual lifetimes of shorter 
than one orbital period close to the last stable orbit. In this case, 
the light-curve will be modulated quasi-periodically at a rate 
close to the frequency of the last stable orbit (see discussion in 
[Eckart et al.|2008a[ and [Eckart et al.|2006b| ). If the hot spots ex- 
pand within the accretion disk or are the source of a short outflow 
above the disk, this quasi-periodic signal will be smeared out. In 
particular, it may be undetectable at lower (sub-)mm observing 
frequencies (see below). In the case of our experiment, only the 
NIR data on 15 May shows a strong sub-flare modulation (as ex- 
plained below) that may be linked to quasi-periodicity. All other 
NIR flares (with the possible exception of the steeply rising and 
falling flanks of the L'-band flare on 26 May 2008) cannot be as- 
sociated with this phenomenon. However, in the case of a source 
expansion this model can be linked to the adiabatic expansion 
model for synchrotron sources that is described in the following 
section. 



3.2. Adiabatically expanding source components 

Adiabatic expansion of the synchrotron components can explain 
the apparent time difl'erence of 1.5+0.5 hours between the NIR 
and sub-mm/ mm flares ( [Eckart et al.|2006a[[2Q08bl[Yusef-Zadeh[ 



et al. 



20Q6a| . 



We assume a uniformly expanding blob of rela- 
tivistic electrons with a power-law energy spectrum, n(E) oc E~p, 
threaded by a magnetic field that declines as R~^, and of energy 
and density that decline as and R~^, respectively, as a result 



of expansion of the blob ( [van der Laan|1966"| ). If Rq, Sq and tq 
are the size, flux density, and optical depth of the source compo- 
nent at the peak frequency (vq) of the synchrotron spectrum, the 
optical depth and flux density at a given frequency v scale as 



Ty = To 



and 



-(p+4)/2 



{- 



-(2p+3) 



5/2 



R\ 1 - exp(-Tv) 
Roj l-exp(-To)' 



(8) 



(9) 



8 
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Fig. 6. Simulated images of multiple spots revolving around a black hole, projected on the observer image plane (a, J3), viewed at 
three inclination angles i=0^, 30^, and 60^, with respect to the common rotation axis. The dotted and dashed circles represent the 
event horizon and the marginally stable orbit, respectively. 




50 lOD tso 

Time [mih] 



Fig. 7. Light curve obtained from the multiple spot simulations shown in Fig. [6j the solid black line representing the overall light 
curve and the different coloured dashed lines showing the contributions of individual spots. 



respectively. Here tq is defined to be the optical depth corre- 
sponding to the frequency at which the flux density is maxi- 
mum, in accordance with van der Laan's model (1966), to be 
able to combine the SSC formalism with an adiabatic expansion 
model.This makes tq dependent only on p through the condition 



e^' - To(p + 4)/5 -1=0. 



(10) 



For instance, if p ranges from 1 to 3, tq ranges from to 0.65. 
Thus for a given particle spectral index p and peak flux den- 
sity 5*0 at a frequency vq, the model gives the variation in flux 
density at any frequency as a function of the expansion factor 
(R/Rq). To convert the dependence on radius to a dependence 
on time, we assume a simple linear expansion model with a con- 
stant expansion speed of Vexp, such that 7?-7?o=Vexp(^^o)• At times 
t<to, we assume that the source has an optical depth equal to its 
frequency-dependent initial value tq at R=Ro. For the optically 
thin part of the spectrum, the flux increases initially with increas- 
ing source size at constant optical depth tq and then decreases 
with decreasing optical depth as it expands. One Schwarzschild 
radius corresponds to 7?^=2GM/c^~10^^ m for a -4x10^ M© su- 
permassive black hole, which infers the velocity of light to be 
about lOORs per hour. For t>to, increasing the turnover fre- 
quency vo or initial source size Rq has the efl'ect of shifting the 



decaying flank of the curve to later times, as does decreasing the 
spectral index a^synch or peak flux density Sq. Increasing the adi- 
abatic expansion velocity Vexp, on the other hand, shifts the peak 
of the light curve to earlier times. Flare timescales are longer at 
lower frequencies and have a slower decay rate, as a result of 
adiabatic expansion. 



3.3. Adiabatic expansion modelling 

Table |3] summarises the properties of the model, such as the adi- 
abatic expansion velocity Vexp, the optically thin spectral index 
o'synch, cutoflf frequency ymax,obs, and flux density 5' niax,obs of the 
source compo nents, using the same nomenclature as |Eckart et al. 
( 2006a 2009| ). The optically thin NIR flux density is represented 
by 5NiR,synch, whilc S NiR,ssc and S x-ray,ssc givc the upper limits 
to the flux densities of upscattered SSC components. A global 
variation in a single parameter by the value listed in the corre- 
sponding column in Table [3] results in an increase of A;^ = 1 
(from the reduced value). Here global variation means adding 
to a single model parameter for all source components the Icr 
uncertainty, such that a maximum positive or negative flux den- 
sity deviation is reached. Alternatively, a variation by the listed 
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date 


model 


source 


^exp 


S max,obs 


Q^synch 


^0 


^max,obs 


B 


S NIR, synch 


S NIR,SSC 


S X-ray,SSC 




label 


hours 


in c 


[Jy] 




[R.] 


[GHz] 


[G] 


[mJy] 


[mJy] 


[nJy] 


1 cr^ 




+ 1.0 


+0.001 


±0.1 


±0.1 


±0.1 


+250 


+ 10 


+ 1.0 


+ 1.0 


+20 



15 May 2007 


A 


a 


0.0 


0.007 


0.6 


0.85 


1.3 


340 


33 


3.1 


<0.03 


<10 






P 


0.0 




0.3 


0.60 


0.3 


840 


29 


12 


<0.03 


<10 






r 


0.30 




0.3 


0.75 


0.2 


840 


32 


5.6 


<0.03 


<10 


15 May 2007 


B 


a 


0.0 


0.017 


0.6 


0.65 


1.3 


340 


37 


8.3 


<0.02 


<40 






P 


0.0 




0.3 


0.65 


0.3 


840 


30 


8.4 


<0.02 


<40 






y 


0.37 




0.3 


0.65 


0.35 


810 


32 


8.4 


<0.02 


<40 


17 May 2007 


A 


a 


0.0 


0.010 


1.0 


1.00 


0.8 


840 


66 


4.7 


<0.01 


<10 






P 


0.45 




0.5 


0.97 


0.3 


1250 


66 


4.0 


<0.01 


15 






y 


-0.59 




1.3 


1.11 


1.0 


840 


69 


3.2 


<0.01 


<10 






6 


-0.96 




0.5 


1.02 


0.3 


1250 


68 


3.2 


<0.01 


12 


17 May 2007 


B 


a 


0.0 


0.011 


0.25 


0.70 


0.2 


1360 


86 


8.1 


<0.01 


19 






P 


0.36 




0.25 


0.70 


0.2 


1360 


86 


8.1 


<0.01 


19 






y 


-0.69 




0.17 


0.70 


0.2 


1360 


65 


5.1 


<0.01 


21 






6 


-1.05 




0.17 


0.70 


0.2 


1360 


65 


5.1 


<0.01 


21 






e 


-0.33 




1.75 


1.05 


1.6 


570 


68 


7.1 


<0.01 


<2 


19 May 2007 


A 


a 


0.0 


0.007 


1.30 


1.07 


1.0 


1360 


67 


3.5 


<0.01 


<10 






P 


2.50 




1.10 


1.10 


0.8 


1360 


67 


0.5 


<0.01 


<10 






y 


2.90 




1.10 


1.10 


0.8 


1360 


67 


0.5 


<0.01 


<10 


19 May 2007 


B 


a 


0.0 


0.017 


1.33 


1.13 


1.0 


720 


68 


2.5 


<0.01 


<10 






P 


0.50 




1.33 


1.09 


1.0 


720 


67 


3.1 


<0.01 


<10 






y 


2.85 




1.50 


1.30 


0.8 


820 


72 


1.2 


<0.02 


<10 






6 


3.30 




1.50 


1.30 


0.8 


850 


72 


1.2 


<0.02 


<10 


26 May 2008 


A 


a 


0.0 


0.005 


1.1 


0.68 


0.5 


1090 


30 


31 


<0.02 


41 






n 
P 


rv /I c 




1 n 
i.U 


U.DO 


U.J 


1 rvnrv 


/in 


ZD 


<U.UZ 


1 o 






y 


0.52 




1.0 


0.83 


0.5 


1030 


31 


15 


<0.02 


17 






6 


1.03 




1.0 


0.80 


0.5 


1160 


34 


17 


<0.03 


23 






e 


1.48 




0.8 


0.77 


0.5 


1160 


48 


16 


<0.01 


<10 








1.90 




0.8 


0.77 


0.5 


1160 


48 


16 


<0.01 


<10 


26 May 2008 


B 


a 


0.0 


0.007 


1.4 


0.70 


0.6 


1090 


44 


39 


<0.02 


24 






P 


-0.40 




1.3 


0.70 


0.6 


1090 


56 


35 


<0.02 


13 






y 


0.52 




1.1 


0.77 


0.6 


1030 


36 


23 


<0.02 


18 






6 


1.03 




1.3 


0.80 


0.5 


1160 


33 


23 


<0.04 


33 






e 


1.48 




1.0 


0.77 


0.5 


1160 


58 


20 


<0.01 


<10 








1.90 




1.0 


0.73 


0.5 


1160 


56 


24 


<0.01 


<10 



Table 3. Source component parameters for the combined SSC and adiabatic expansion model of the 15, 17, and 19 May 2007 and 
the 26 May 2008 flares. Labels A and B refer to models with lower and higher expansion velocities, respectively. The flare times 
are given with respect to the peak of the brighter NIR flares. In addition to Vexp, the 7?o values are responsible for the position and 
width of the infrared flare peaks in time. Different values for Ofsynch are required to match the infrared flux densities. 



uncertainty for only a single source component results in a varia- 
tion of the model predicted NIR and X-ray flux density by more 
than 30%. Judging from the l^x based on the mm-data alone, the 
global uncertainties for 5'niax obs, Q^synch, and may then be dou- 
bled. 

When performing the reduced fit for N source components, 
we used times 4 (5'niax,obs, o^synch, ^o, Vmax) plus one com- 
mon expansion velocity Vexp and time ofl'set (leaving the time 
differences between the components fixed), i.e., AN +2 degrees 
of freedom. Unfortunately, the model parameters may not all be 
considered independent, e.g., the width and peak of a light curve 
signature depends to a varying extent on all 4 parameters 
<^synch, ^0, and Vmax- Therefore, for all models we stayed with the 



minimum number of source components to estimate the degrees 
of freedom. Since in the case of the May 2007 flare the VLT 
data consist of -10 times the number of CARMA data points, 
we weighted the squared CARMA flux deviations and number 
of data points by an additional factor of 19. The;^^ test was then 
carried out using the sum of the squared flux deviations and data 
points of the VLT and CARMA datasets. 
In general, we aimed to develop a model at a lower and a higher 
expansion speed labelled A and B (Table [3]). In comparison to 
earli er modelli ng result s (Eckart et al. 2QQ8a , 2009 ; Yusef-Zadeh] 
|et al. 2008 ; M arrone et al. 2008 ), these speeds were chosen to be 
0.007c and 0.017c for the 15 May and 19 May 2007 flares. For 
the 17 May 2007 and the 26 May 2008 data, a violation of the ob- 
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served flare width and the expected range for the magnetic field 
strength restricted the choice of expansion speeds. For compar- 
ison, we listed the model components in Table [3] The models 
were developed by minimizing the number of free parameters 
(and maximizing the description of significant flare features in 
the observed light curves). 

SSC modelling as an additional constraint: The SSC model 
described in Sect. 13.11 allows us to estimate the SSC contribu- 
tion of the X-ray and NIR flux densities, and the magnetic field. 
Higher SSC NIR emission would violate the assumption in our 
adiabatic expansion model that the dominant source of NIR flux 
density is synchrotron emission of the THz peaked expanding 
component, as well as require us to model the un known X-ray 
flux density. Magnetic fields of up to -70 gauss ( Eckart et al. 
|2QQ6a| [2Qn8bl |Yusef-Zadeh et al.||20Q8l [Marrone et al.||2QQ8| ) 
have been obtained for previous models via B ~ G^v^^S'^, where 
Vm and S;;^ are the synchrotron turnover frequency and flux den- 
sity, respectively. These results should also not be violated by 
the model. 

The expansion speed: The model yields expansion velocities 
from 0.005c - 0.017c (see Table [3|, which agree well with pre- 
viously published results ( Yus ef-Zadeh et al.||2008| . These ve- 
locities are lower than expected for relativistic sound speed in 
orbital velocity near the SMBH. This may be caused by the bulk 
velocity of the source components being larger than the expan- 
sion velocity Vexp or because the expanding material is confined 
to the immediate region surrounding Sgr A* in the form of a 
disk or corona. In this case, shearing caused by diff'erential rota- 
tion within the accretion disk is responsible for the 'expansion' 
and the observed low expansion velocities ( Eckart et al.|2008b{ 
Zamaninasab et al.|2008bl|Pechacek et al.|2 008 ). Owing to the 



expansion or the presence of several strong spots in the disk or 
corona, the short time modulation often seen in the NIR will also 
be smeared out in the (sub-)mm- wavelength domain - implying 
that it may not be observable at all at these wavelengths. 
By comparing other flare events that were detected simultane- 
ously at X-ray and sub-mm/mm-wavelengths, we imposed the 
additional constraint that the 100 GHz flux density should not 
be much stronger than about 0.2 Jy, 1.5+0.5 hours after the NIR 
event, and the X-ray flux density should be limited to a few 
lOnJy. The following sections describe the modelling of the in- 
dividual flares. 

3.3.1 . Expansion modelling of the 1 5 May flare emission 

The 15 May 2007 flare can only be constrained by means of its 
NIR lightcurve (Fig.jSj) since the mm data begins nearly 3 hours 
after the first peak in the NIR. According to the adiabatic expan- 
sion model, 3-4 hours after the NIR flare, the flux density in the 
mm band was insuflSciently strong to be significant compared 
to the median flux density variations in the data. We modelled 
a background component p and two components for the two 
strongest sub-flares a and y. The expansion speed could not be 
constrained at all. In model A, we used a value of 0.007c, close to 
the speed we found for most of the modelled flares. In model B, 
we used a higher speed of 0.017c, which required lower spectral 
indices and a weaker background component to obtain a good fit 
to the NIR data and to fullfil the upper limits at 100 GHz and in 
the X-ray domain. 

A correlation between the L'-band data at 10 UT and the ra- 
dio data cannot be excluded, but the L'-band data light curve is 
not long enough to justify including this data in a full modelling 
of the NIR/mm-light curves in the framework of the adiabatic 
expansion model. However, high NIR flux density levels pre- 
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Fig. 8. Models A and B of the 15 May 2007 NIR flare: a) Fitted 
NIR light curve, with the black curve corresponding to the model 
and the red curve corresponding to the data, b) Contribution 
of individual source components to NIR light curve, and c) 
Contribution of individual source components to mm light curve. 



cede a positive flux density excursion in the mm-domain, as it is 
expected to in the framework of the adiabatic expansion model. 
Although the L'-band flux density is higher than that at other 
times, the event is not linked to a particularly strong flux den- 
sity excursion in the mm-domain. However, the high L'-band 
flux may also be caused by a transient steeper infrared spectral 
index. This may indicate either a significantly higher turnover 
frequency (extrapolating into the sub-mm domain) or the pres- 
ence of stronger synchrotron losses. Both eff'ects may explain a 
lack of strong mm-emission following the NIR event at 10 UT. 

3.3.2. Expansion modelling of the 17 May flare emission 

The 17 May 2007 flare is constrained by the NIR and CARMA 
mm-data (Fig. [9]). This time lag of 1.5+0.5 hours between the 
flare events mainly determines the expansion velocity. A con- 
stant flux density of 0.003 Jy was subtracted from the NIR data 
to model the individual flare events. We present two modelling 
approaches that lead to a successful description of the flare. In 
model A, we used only 4 components describing 4 sub-flares a 
through d. Motivated by the relativistic disk modelling referred 
to in Sect. |3.1| we used in model B a background component 
6 located close to the centre of the overall NIR flare. We then 
modelled the 4 sub-flares with additional components. In com- 
parison to model A, this results in lower values for 5'niax,obs and 
the spectral index CKsynch, smaller source component sizes 7?o, and 
higher cutofl' frequencies Vmax,obs for components a through 6. 
The quality of the fit is comparable for models A and B. 
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Fig. 9. Models A and B of the 17 May 2007 NIR flare: a) Fitted 
NIR light curve, with the black curve corresponding to the model 
and the red curve corresponding to the data, b) Contribution of 
individual source components to NIR light curve, c) Fitted mm 
light curve, with the black curve corresponding to the model and 
the red curve corresponding to the data, and d) Contribution of 
individual source components to mm light curve. 



Fig. 10. Models A and B of the 19 May 2007 NIR flare: a) Fitted 
NIR light curve, with the black curve corresponding to the model 
and the red curve corresponding to the data, b) Contribution of 
individual source components to NIR light curve, c) Fitted mm 
light curve, with the black curve corresponding to the model and 
the red curve corresponding to the data, and d) Contribution of 
individual source components to mm light curve. 



3.3.3. Expansion modelling of the 19 May flare emission 

For the 19 May 2007 flare, the slower adiabatic expansion ve- 
locity of 0.007c allows us to fit the first event (centred on the 
NIR at -5:30 UT) with a single component (Fig. 10). For higher 
velocities, a larger number of source components has to be used 
to model the NIR data. For both velocities, the second, later flare 
component (centred on the NIR at -8:20 UT) is constrained 
by the 3 mm CARMA data. For the model parameters used to 
match the NIR data, the predicted 3 mm flux densities of both 
flare components (at -5:30 UT and -8:20 UT) are at or below 
the flux density limit of 0.1 -0.2 Jy, above which no millimeter 
wavelength flare event has been detected. 



3.3.4. Expansion modelling of the 26 May flare emission 

Modelling the 26 May 2008 flare (Fig. [TT]) requires a satisfac- 
tory description of an initially bright flare with sharply rising and 
falling flanks followed by a 1.5 hour plateau. This can only be 
achieved with several source components each covering a max- 
imum portion of time of about a 30 minute duration. A smaller 
number of source components would require larger source com- 
ponent sizes resulting in higher magnetic fields and a violation 
of the 100 GHz flux limit. Given these difficulties, the number 
of source components is of course not unique and the models 
proposed in Table [3]can only be taken as an example. 

The overall flare structure is quite reminiscent of the flare 
that was observed on 3 June 2008 simultaneously using APEX 
and the VLT, as described by Eckart et al. ( 2008b ), with the ex- 
ception that the 26 May flare is more continuous than the 3 June 
2008 flare. An overall flare structure that starts out with a bright 
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Fig. 11. Models A and B of the 26 May 2008 NIR flare: a) 
Fitted NIR light curve, with the black curve corresponding 
to the model and the red curve corresponding to the data, b) 
Contribution of difl'erent source components to NIR light curve, 
and c) Contribution of diff'erent source components to mm light 
curve. 



event foll owed by a single or a few fainter events was also re- 
ported by lEcka rt et al. | ( 2006a| 2009^). These events may relate 
to a common physical structure. An explanation could be a disk 
structure that is expanded by difl'erential rotation within the disk. 
Alternative explanations may involve a common history in terms 
of either accretion or m agnetic field instability - possibly also 
within a jet structure (see [Eckart et al.| ( [2008b |a| ) and Sect.|4]). 



4. Alternative short jet model 

Emission from a jet and an underlying accretion process may 
also account for the spectrum of Sgr A*, with the plasma in the 
jet starting to become optically thin at longer wavelengths with 
increasing distance from the centre. Thus, radiation at difl'erent 
radio wavelengths probe difl'erent sections of the jet, leading to 
a correlation with emission in other wavelength regimes. Very 
close to the SMBH, it is difficult to distinguish between emis- 
sion from the jet and the accretion flow. For this jet model, flux 
variations may then be caused by accretion processes or jet in- 
stabilities rather than be the result of modulation from an orbit- 
ing spot, followed by an adiabatic expansion of jet components. 
Details of a comp act, weak jet st ructure are discussed in [Markoff 
[etaL] ( [20071 ) (also [Markoff eTa". tlOOSj . 



5. Summary 

We have presented results from global coordinated multiwave- 
length observational campaigns carried out in 2007 and 2008 
using NACO at VLT in the NIR K- and L' -bands, and CARMA, 
ATCA, and IRAM in the mm regime. 

We have presented a new method to obtain concatenated 
light curves of the compact mm-source Sgr A* from single dish 
telescopes and interferometers in the presence of significant flux 
density contributions from an extended and only partially re- 
solved source. The method requires several consecutive datasets 
to be available for each participating observatory. From these 
data, we have evidence of four different ff aring events in the NIR, 
with three of these events being covered later in the mm regime. 

We have modelled the ffare emission with a model involving 
adiabatic expansion of synchrotron source components, and ob- 
tained spectral index values, expansion velocities, and other pa- 
rameters that are consistent with previously published variabil- 
ity studies. The mm-ffares were found to occur on average about 
1.5+0.5 hours after the NIR ffare emission, and the expansion 
velocities for the various ffares ranged from 0.005c - 0.017c. 

Modelling the NIR-ffares with synchrotron components that 
become optically thick in the sub-mm wavelength regime and 
obey the ffux density values or limits in the mm- and X-ray do- 
main constrains their turnover frequency and ffux density. For 
the described global experiment, the low mm-regime ffux den- 
sity limits imply that the turnover frequencies are not lower 
than about 1.0+0.3 THz and that the optically thick peak ffux 
densities at or below this turnover frequencies do not exceed 
about ~lJy on average. However, in other experiments the 
source showed higher ffux density ffare emission in the mm- 
domain (e.g., 0.5-1.5 Jy in the sub-mm; Yusef-Zadeh et al. 2009, 
Marrone et al. 2008, Eckart et al. 2008, 2006b) 

Further monitoring of the variability of Sgr A* in diff'erent 
wavelength regimes (including X-ray, NIR, sub-mm, and mm 
radio) and in polarised NIR/radio emission is required to im- 
prove our understanding of the adiabatic expansion model and 
improve statistics describing the ffaring activity. Current mm- 
interferometers such as CARMA, PdBI, and ATCA and future 
telescopes such as ALMA which are able to distinguish between 
emission from Sgr A* and the thermal emission of the CND and 
mini-spiral surrounding the SMBH, and NIR telescopes with 
large apertures such as VLT, Keck, and LBT, which can dis- 
tinguish SgrA* from surrounding stars, may be combined to 
provide us with high quality data to study the evolution of syn- 
chrotron components in greater detail. With future mm- VLB I 
measurements at frequencies of 230 GHz and above, imaging of 
the central region of Sgr A* may become possible, enabling a 
deeper study of the accretion p hysics and testing of current the- 
ories of emission in the region ( Doeleman et al.|2008[ 2009| ). 
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Appendix:NIR light curves 

The individual NIR K and L'-band light curves from Fig. [T] are 



shown in Figs. 12 to 16 A spectral index of -0.9 (as we preferen- 
tially used in our modelling, see Table[3]and Eckart et al. 2008b, 
2009, Yusef-Zadeh et al. 2008) gives us a scaling factor of 0.6, 
which we use to scale the L'-band data, and then combine the K 
and L'-band data as shown in Fig.[T]and Figs.[T2]l 



to 



16 In the re- 



maining figures, no scaling has been applied. The uncertainties 
in flux measurements were obtained from the reference star S2, 
which has known flux densities of 22+1 mJy and 9+1 mJy in the 
K and L' -bands, respectively. 
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Fig. 16. Observing time-K-band: 04:55:00 to 09:28:22 
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